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^ ; Abstract 

A theoretical study is presented to quantitatively analyze the transport ex- 
periment through individual DNA molecules reported recently by Porath et 
al. [Nature 403, 635 (2000)]. A variety of valuable quantities are identified 
q ■ by contacting the theoretical model with the measured data. The partially 

o ; 

decoherent nature on the GC pairs of DNA is elaborated in contrast to the 
^- ■ completely incoherent hopping mechanism discussed in the context of charge 

transfer experiments. 
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DNA has the special double-helix structure with 7r-electron cores of well stacking bases, 
which may provide an one-dimensional pathway for charge transport. Since first posed in 
1960s, the question whether or not DNA is able to conduct electrical charges is still debated. 
Due to the advent of molecular electronics as well as other biological considerations, this issue 
has stimulated intense research interest. Recent studies proposed that DNA seems to be one 
of the most promising candidates in the application of functional nano-electronic devices [|TJ . 
Very recently, direct measurements for the electrical transport through DNA molecules were 
performed 10,0]. In particular, Porath et al. |§ measured the electrical transport through 
the individual DNA molecules. The observed nonlinear current- volt age (I-V) characteristics 
clearly suggest that DNA molecules be good molecular semiconductors. In this letter, we 
carry out a quantitative analysis for the experimental I-V characteristics ||, and discuss the 
relevant physical implications. 

The theoretical model for the experimental setup is schematically shown in Fig. 1. In the 
experiment ||, the DNA sample of 10.4 nm-long chain containing 30 GC pairs was studied. 
The ab initio calculation || showed a large HOMO-LUMO gap (10.07 eV) of the single 
GC pair. We thus adopt a homogeneous one-band tight-binding model for the HOMO- 
mediated charge transport through the DNA molecules under study. To proceed, consider 
the following Hamiltonian to describe the system in Fig. 1: 

H = H + Hi, 
H = H M + H L + H R + H res . (1) 

Here H describes the separated subsystems of the DNA molecule (Hm), the left and right 
metal electrodes (Hl and Hr), and the dephasing reservoirs (H ieB ). Hi couples the DNA 
molecule to the electrodes and the dephasing reservoirs. Following Biittiker's idea for the 
inelastic scattering effect [|,0, the phase-breaking processes on the GC-pairs are modeled by 
coupling each GC-pair to & fictitious electronic reservoir, see Fig. 1. The dephasing strength 
is characterized by the coupling parameter 77. 

To carry out the nonlinear transport characteristics, the knowledge of the energy- 



dependent transmission coefficient through the molecule is required. Owing to the open 
boundary conditions associated with the voltage electrodes and dephasing reservoirs, the 
Green's function method would be the proper and versatile theoretical tool. We thus intro- 
duce the following Green's function 

g(E) = (E-n cS )-\ (2) 

where Ti e ^ is the effective Hamiltonian of the reduced molecular system, obtained by elim- 
inating the degrees of freedom of the two electrodes and the dephasing reservoirs. As a 
simple model, the electrodes and dephasing reservoirs can be described by semi-infinite 
one-dimensional tight-binding chains ||, which result in the self-energy corrections to the 
Hamiltonian of the molecule as 

N 

H eS = H M + Z L \l){l\+Z R \N){N\ + '£Z j \j){j\. (3) 

i=i 

El(r) and Ej are respectively the self-energies resulting from the coupling of the molecule 

to the left (right) electrodes, and the j'th dephasing reservoir. Using Dyson equation, the 

self-energy E^ (/i = L, R, and{j}) can be carried out as 

where a^ is the self-energy correction of the semi-infinite chain to the ending-site of the 
chain attached to the DNA molecule, and has the following expression H 
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Here V^ are the coupling strengths between the molecule and the electrodes/reservoirs: 
Vfj, — Vl(r) if /i = L(R); and V^ — r\ if n — {j}. E^ and 7^ denote the site energies 
and the nearest-neighbor hopping strengths of the homogeneous tight-binding electrode and 
reservoir chains. We let 7l(h) = 7 for the left (right) electrode and 7-,- = 70 for the dephasing 
reservoirs. Thus, 47 and 470 characterize the band widths of the electrodes and the dephasing 
reservoirs, which are assumed in this work to be commonly 5 eV. 



To connect the Green's function with the transmission coefficients between any pair of 
reservoirs, say, from fi to v, consider a linear transport between them. It can be shown |7j 
that the linear conductance at low temperatures is g^ = (2e 2 / nh) A,j,A u \Q ^(E)] 2 . Here, A M 
is associated with the self-energy by A M = — ImS M . Q lxv {E) is the matrix element of G(E) 
between the molecular sites connecting with the /xth and vtYi reservoirs (//, v = L, R, andjj}). 
On the other hand, the linear conductance is given by the Landauer-Biittiker formula || 
as g^u = (e 2 /h)T lu/ . We therefore establish the following desired relationship between the 
transmission coefficient T^ u and the Green's function: 

7V = 4A A A,|£v(£)| 2 . (6) 

Note that in the wide band limit, A^ ~ V 2 / r y fl . We thus conclude that in the special case 
Vp — 7/i j Eq. @ reduces to the result of Refs. [§,§. After knowing T^ u , we are ready to 
carry out the total effective transmission coefficient from the left electrode to the right one. 
We have recently shown that for each individual electron being emitted from the electrode, 
the effective transmission coefficient in the presence of phase-breaking is of the same form 
as that first derived by D'Amato and Pastawski JiT]|,p[]: 

N 

% S (E) =T lr +J: kWw^kW (7) 

Here K^f' = T^, KJ) R ' = T vR . W~ x is the inverse of the matrix W defined in terms of 
the matrix elements W^ = (1 - -R W )<V ~ ^(1 - <W> witn R ^ = * _ ^W") T vli . The 
physical meaning of Eq. (|7|) is clear: the first term is the coherent contribution of tunneling, 
whereas the second term denotes the incoherent component owing to electron suffering the 
dephasing processes. 

With the knowledge of the transmission function %s(E), it is straightforward to evaluate 
the I-V characteristics by applying the standard formalism based on the scattering theory 
of transport || 

9p r°° 
I = ^J_JE% n {E)[f L {E) - f R (E)}. (8) 



Here f L / R (E) = {exp[(E — //£/#) /fceT] + 1} _1 is the Fermi function. The room temperature 
(fc B T ~ 26 meV) will be considered in the numerical results, in order to keep consistency 
with the experiment [[J. jjll and fiR are the electrochemical potentials of the two metal 
electrodes, whose values depend on the applied bias voltage. Following Datta et al. ||11|| , 
we set //£, = Ef + (1 — ft)eV, and hr = Ef — neV, where Ef and V are respectively the 
equilibrium Fermi energy and the applied voltage, k is a useful parameter in characterizing 
how the applied voltage V is divided across the two junctions between the molecule and the 
electrodes. 

Figure 2 shows the calculated I-V characteristics, together with the experimental data 
taken from Ref. ||. Below we discuss how to make contact of the theoretical model with 
the measured data in order to get useful informations. Firstly, we notice the asymmetric 
gaps V}^ under positive and negative bias voltages: V^ 1 = 1.6 eV and V}~^ = 0.8 eV. 
This asymmetry suggests k = 1/3, instead of 1/2 as usual [|TTJ. Secondly, the experimental 



differential conductance || showed a peak structure with the voltage spacings AV of 0.1 ~ 
0.5 eV. By taking the average value AV ~ 0.25 eV and k = 1/3, we estimate the DNA 
molecular energy level spacing as AE = kAV ~ 0.08 eV. Since the DNA under study 
consists of 30 GC-pairs, from the simple consideration that the energy bandwidth 4t = 
30 x 0.08 = 2.4 eV, we obtain the important information for the hoping strength between 
the nearest-neighbor GC pairs, t ~ 0.6 eV, which is in good agreement with the ab initio 
calculation JJ. Thirdly, the equilibrium Fermi level Ep of the electrodes can be estimated 
via E F -E G = 2t+KV c (+) ~ 1.73 eV, where E G is the HOMO level of the single GC pair. This 
information would be helpful in understanding and predicting the voltage gaps in transport 
experiments. Finally, we notice the asymmetry of the entire profile of the experimental I-V 
curve under positive and negative bias voltages. We temporarily understand this feature by 
the distinct coupling strengths of the DNA molecule to the metal electrodes: Vl/r = 3 meV 
for positive voltage, and Vl/r = 2 meV for negative voltage. With these parameters, the 
theoretical results are in quantitative agreement with the measured data as shown in Figure 
2. 



Since the I-V characteristic is an integrated result of the transmission coefficient, the de- 
phasing effect is not quite prominent on the I-V curve, except that the voltage gaps and the 
current steps due to the discrete molecular energy levels are smoothed by the phase-breaking 
scattering as shown in Fig. 2. However, it has drastic effect on the differential conductance 
as shown in Fig. 3. Roughly speaking, the dephasing would broaden the conductance peaks. 
More precisely, the differential conductance dl/dV which corresponds to the transmission 
coefficient T e s(E) at the specific voltage would change with the dephasing strength as rep- 
resentatively shown in the inset of Fig. 3. The turnover behavior has also been discovered in 



other context of tunneling with dissipation ||10| . Note that the experimental peak width is 



of ~ 0.2 eV, while its spacing ~ 0.25 eV. It may thus suggest that the charge motion along 
the stacked GC-pairs is only partially dephased in the referred experiment [[|. It is worth 
emphasizing that this motional nature is essentially different from the hopping mechanism 
discussed in the context of long-range charge transfer phenomena in DNA [0], where the 



quantum coherence is regarded as being completely destroyed on the GC pairs. We may 
attribute this distinct nature to the different environments in which the DNA molecules 
located. Since the electronic dephasing would significantly influence the long-range charge 
transfer efficiency [|IUJ, we here suggest to further identify the dephasing strength on the GC 
pairs, by investigating the electrical transport through various sequence of DNA molecules 
and under different environments. 

In summary, we have presented a theoretical formalism for the transport through DNA 
molecules. By making contact of the theoretical model with the experimental data, a va- 
riety of valuable quantities were identified, and the theoretical result was in quantitative 
agreement with experiment. The modest dephasing nature on the GC pairs was emphasized 
in particular, viewing its close relevance to the long-range charge transfer in DNA. Finally, 
the present theoretical formalism is quite general, which can also be employed to describe 
the transport through other molecular wires. 
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FIGURES 
FIG. 1. Schematic diagram for the electrical transport through DNA molecules. The DNA 

molecular wire consists of a stack of GC pairs, and each GC pair is connected with a fictitious 

electronic dephasing reservoir by a coupling strength r\. 

FIG. 2. I-V characteristics: theoretical results versus experimental data taken from Ref. [3]. 
The transport currents in the presence of weak dephasing (77 = 0.05 eV) and stronger one (77 = 0.3 
eV) are theoretically shown by the solid and dashesd curves. 

FIG. 3. Dephasing effect on the differential conductance. The solid and dashed curves corre- 
spond to 77 = 0.05 eV and 7/ = 0.3 eV, respectively. The detailed change of the conductance profile 
is representatively shown in the inset. 



Fig. 1 by X.Q. Li et al: Electrical transport . . . 
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Fig. 2 by X.Q. Li et al: Electrical transport 
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Fig. 3 by X.Q. Li et al: Electrical transport 
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